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ABSTRACT 

The shaping of the nebula is currently one of the outstanding unsolved prob- 
lems in planetary nebula (PN) research. Several mechanisms have been proposed, 
most of which require a binary companion. However, direct evidence for a binary 
companion is lacking in most PNs. We have addressed this problem by obtaining 
precise radial velocities of seven bright proto- planetary nebulae (PPNs), objects 
in transition from the asymptotic giant branch to the PN phases of stellar evolu- 
tion. These have F-G spectral types and have the advantage over PNs of having 
more and sharper spectral lines, leading to better precision. Our observations 
were made in two observing intervals, 1991-1995 and 2007-2010, and we have 
included in our analysis some additional published and unpublished data. Only 
one of the PPNs, IRAS 22272-1-5435, shows a long-term variation that might ten- 
tatively be attributed to a binary companion, with P > 22 years, and from this, 
limiting binary parameters are calculated. Selection effects are also discussed. 
These results set significant restrictions on the range of possible physical and 
orbital properties of any binary companions: they have periods greater than 25 
years or masses of brown dwarfs or super- Jupiters. While not ruling out the 
binary hypothesis, it seems fair to say that these results do not support it. 

Subject headings: binaries: general - binaries: spectroscopic - circumstellar matter - 
planetary nebulae: general - Stars: AGB and post-AGB - Stars: mass-loss 
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1. INTRODUCTION 

1.1. Background to the Question 

Arguably the most controversial area in planetary nebula (PN) research at present is 
the determination of the mechanism for shaping the nebula. This has been brought to the 
fore by the visually stunning, high-resolution Hubble Space Telescope (HST) images. PNs 
generally possess an ellip tical or bipolar struct ure, often with additional point symmetric 



features (ansae, jets; see 



Balick fc Frank 



2002). In contrast, their precurs ors, asymptotic 



giant branch (A GE) stars, have been considered to be basically spherical ( lOlofsson 



Neri et al. 



1999 



19981 ). although some recent resolve d molecular-line observatio ns show that a 



20071 ). So the question 



significant fraction show some axial symmetry (iCastro-Carrizo et al. 
can be posed as to how a mass-losing AGE star, which is basically spherical, can evolve 
into the variety of PN shapes. 

The presently emerging consensus is that the axial and bipolar asymmetry is caused 
directly or indirectly by the presence of a binary companion to the central star. A binary 
companion could influence the density structure in at least the following three ways. 



1. Most directly, a companion could gravitationally focus the mass loss into the orbital 
plane, forming an equatorial density enhancement and perhaps a torus. This density 
enhancement or toru s would then collimate the fast wind producing a bipolar outflow 



(ILivio fc Soker 



19881 ). A variation on this would be the formation of an accretion disk 



around the binary companion; this could collimate a fast wind and carve out the lobes 



and also lead to point-symmetric ejecta (ISoker &: Rappaport 



2000). 



2. The mass could be preferentially lost in the equatorial plane of the PN during the 



AGE phase due to a rotational. 



torus (iGarcia-Segura et al. 



y-induced oblate shape, producing a coUimating 



19991 ). However, given the low rotational velocity of an 
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intermediate-mass star on the main sequence, it would probably require the presence 
of a close companion to spin up the resulting AGB star by the transfer of orbital 
angular moment um. This binary-induced ma ss loss could occur during a common 



envelope phase (INordhaus fc Blackman. 



20061). 



3. The c entral star could possess a 



lobes (iGarcia-Segura et al. 



magne tic field that coUimates the outflow into bipolar 



1999 



20051 ). Recent work indicates that something like a 



binary interaction would most probably be needed to sustain such a magnetic field 



through the transfer of angular momentum (INordhaus et al. 



20071 ). 



The first mechanism is generally favored to produce the density enhancement or 
collimating torus, although it can be seen that a binary companion would be important 



in all three mechanisms. Population synthesis stu dies suggest that 



jinary interactions 



20061 ). Thus it is 



can produce the correct number of Galactic PNs (iMoe fc De Marcol 
increasingly common to hear it stated that elliptical or bipolar PNs are due to the effect 
of a binary companion, or even that the presence of an axially-symmetric or bipolar 
structure implies a binary. However, this hypothesis has not been adequately supported 
observationally. An extensive review of the question of what shapes PNs, with a thorough 
investigation of the evidence for the binary hypothesis, has recently been presented by 
De Marcol (120091 ): the need for direct observational tests of this hypothesis is stressed 



^To address this outstanding problem in the study of PNs, an informal international 
collaboration was recently formed with the goal of settling the question of the binarity 
in PNs and its effect on the shaping of the nebula (PlaN-B; coordinator O. De Marco; 



http:/ /www. wiyn.org/plan-b/ ). 
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1.2. Why Search for Binary PPNs? 



Proto-planetary nebulae (PPNs) are the objects in transition between the AGB and 
PN phases in the evolution of intermediate- and low-mass stars. During the AGB phase, 
such stars are surrounded by an expanding circumstellar envelope CSE of mass being 
lost at an increasing rate. In the PPN phase, the high rate rnass lo ss has ended and the 



star is surrounded by a detached, expanding envelope (IKwok 



19931 1. PPNs are a subset 



of the larger class o f post-AGB stars, which also includes RV Tau and R CrB variables 



( iVan Winckel 



20031 ). PPNs can be distinguished by having more massive circumstellar 



nebulae and in many cases showing clear abundance patterns from AGB nucleosynthesis. 
As such, PPNs appear to be the most likely objects to evolve into PNs. In the Discussion 
section, we will distinguish in more detail PPNs from post-AGB stars in general and 
compare our results with those found for other post-AGB stars. 

PPNs display a basic axial symmetry, often showing bipolar lobes and occ asionally 
point symmetry. This has been particularly seen in high— res olution i7>ST images ( lUeta et al. 



2000 



Su et al. 



2001 



Sahai et al. 



2007 



Siodmiak et al 



20081 ) . Some also display an obscured 



equatorial region. Thus one sees in PPNs the same basic structures as in PNs, but at an 
earlier stage in the nebula, a stage closer to the beginning of the shaping process. This 
commends the study of PPNs to investigate the shaping mechanism(s). 

Binarity can manifest itself in several ways: a visible companion, photometric light 

velocity variat ions. A survey of the results of 



2009al ) with null results; no evidence of 



variations, composite spectrum, and radia 
these methods has been presented earlier (iHrivnak 
binary companions to PPNs was found. 

In this paper, we discuss the observational evidence of binarity in PPNs based on 
long-term radial velocity studies of seven bright PPNs. The results of this study and their 
implications for the binary nature of PPNs are then discussed, and conclusions drawn and 
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discussed on whether or not they provide evidence to support the binary hypothesis. 



OBSERVATIONAL SAMPLE AND DATA SETS 



Our sample consists of the seven brightest PPNs observable from mid-latitudes in the 
northern hemisphere. They range in visual brightness from 7.1 to 10.4 mag and are listed 
in Table [H We examine later the question of what biases might result from observing the 
brightest PPNs. For uniformity, we will refer to all of them by their IRAS identification 
numb ers. All show the d ouble-peaked spectral energy distribution characteristic of PPNs 



[e.g. 



Hrivnak et al. 



19891 ). and all have F— G spectral types. All are known to vary in light 



and all but one have variable star desig nations. Detailed light curve studies of all seven 



of these hav e been presented elsewhere ( iHrivnak et al 



Fernie 



19831 ) 



20101 : 



Arkhipova et al 



2010 



2006 



EDITOR: PLACE TABLE □ HERE. 



The main data sets used in this study are the radial velocity observations we carried 
out at the Dominion Astrophysical Observatory (DAO), initially from 1991-1995 and then 
more recently from 2007-2010. The 1991 -1995 data were obtained with the radial-velocity 



spectrometer (RVS, 



Fletcher et al. 



19821 ). while those from 2007-2010 were derived from 



CCD spectra by cross-correlation with several lAU standards. However, the velocities of 



those standards w ere taken to be those derived at the DAO photographically ( IScarfe et al. 



1990 



Scarfe 



20101 ) and with the RVS. The zero-point of the much more numerous RVS 



observations of standard stars has been adjusted to match that of the photographic data 
for this purpose, this ensuring that the CCD and RVS data are on the same system. With 
F— G spectral types, these objects have numerous sharp absorption lines, which result in an 
observational precision of ~0.7 km s^^. The individual radial velocity measurements will 
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be published elsewhere in the context of our detailed pulsational studies of the individual 



objects. We find that they a 



ranging from 35 to 130 days (IHrivnak 



] vary in 



velocity and light due to 



(12000 



2009bl : 



Hrivnak et al 



pulsation , with periods 



2010. see also 



Arkhipova et al. 



20101 )). Substantial radial velocity data sets have also been publis. 



inves tigators for thr ee of the ob 



2000), 17436+5003 f,Burki et al. 



ects, IRAS 07134+1005 flLebre et al 



19801), and 22272+5435 fiZacs et al. 



l ed by other 



1996 



Barthes et al 



20091 ). as part of 



their study of pulsation in these objects. These published data sets and a few individual 
published and unpublished velocity observations have also been incorporated into this 
study. 



3. RESULTS OF THE BINARY INVESTIGATION 

3.1. Investigating the Sample 

In Figures [Tal and [Tbl are plotted the radial velocity data for these seven objects. In the 
left-hand panels are shown our observations from 1991-1995 on a scale that allows one to 
get a better sense of the velocity variations, and in some cases the clear cyclical nature of 
the pulsations is indicated by sine-curve fits. In the right-hand panels are shown the data 
over the entire range of observations, including those by others. The combined data sets for 
each of the objects show a range of velocities of 10 to 14 km s~^. 

EDITOR: PLACE FIGURE [li HERE. 



EDITOR: PLACE FIGURE [Tb] HERE. 



Comparing our 1991-1995 data with 2007-2010 data for the seven objects (summarized 
in Table [2]), we find only one object, IRAS 22272+5435, for which the average values differ 
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by more that 1.5 km s~^(or 2a). For the others, the difference between the two data sets 
is < 2a and it does not appear that the average values differ significantly between the two 
epochs of observation. Nor is there evidence for systematic change when we include the 
other data sets. A formal period analysi s of all of the data was carried out for each object 



using the Period04 ( iLenz fc Breger 



l2005l ) program. Beyond the pulsation mentioned earlier, 
there is a suggestion of a long-term periodicity only in IRAS 22272+5435, which will be 
discussed in Sec. 13.21 



EDITOR: PLACE TABLE M HERE. 



To investigate even longer-term variations, we compared the average radial velocity 
from these optical, photospheric spectra with the radial velocity determined from the 
circumstellar CO or OH emission. The CO and OH represent the emission from the 
circumstellar envelopes, which have been expanding over several thousand years, since the 
late stages of the AGB phase. Thus their velocity centers should be the same as that of 
the average PPN velocities, which should be that of the optical PPN velocity if there is 
no binary motion and that of the barycenters if there is binary motion. These results are 
also listed in Table O We see that these molecular-line velocities, when transformed to the 
heliocentric system, are all, with one exception, close to the velocities based on the visible 
spectra, differing by to —2 km s~^. The lone exceptio n is IRAS 18095+ 2704, which differs 



by +5 km s ^. This is the only target detected in OH ( lEder et al. 



19881 ). and inspection of 



the spectrum shows that the two separate maser components are relatively weak compare d 



Bujarrabal et al. 



fll992h 



to most of the others in the study and that they show structure, 
cite a tentative CO detection at a velocity that would reduce the difference to +2 km s~^. 
Based on these data, we do not see any evidence from the molecular-line observations to 
indicate a significant difference between the photospheric and the circumstellar envelope 
velocities, although there appears to be some inconsistency with the OH velocity of IRAS 
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18095+2704. 



3.2. Evidence for a Possible Binary Companion to IRAS 22272+5435 



A comparison of the earlier (1989 -1995) anc 



22272+5435, including the data from iZacs et al 



ater (2005-2010) observations of IRAS 



(120091 ). shows a bimodal distribution of 



velocities. There exists a significant difference of —2.2 km s~^ in the average velocities of 
the data sets from the two epochs, about four times the sum of the uncertainties in the 
values. We interpret this velocity difference as most probably due to the motion of the 
PPN around the barycenter of a binary star system. We note that IRAS 22272+5435 is the 
coolest and reddest star in the sample and presumably has the lowest surfac e gravity, as 



indic ated by it spectral type and the results of a model atmosphere analysis (IReddy et al 



20021 ). Hence it is likely to be the least stable on a time scale of many years. While we have 
considered this possibility, we still think that the binary hypothesis is more likely to explain 
the long-term velocity variation. 

A formal period study was carried out of all of these data. In addition to the pulsation 
period of 128.3 ± 0.1 days, no reliable long-term period could be determined. This is not 
surprising, give the distribution of the velocities in the two observing intervals around two 
nearly constant average velocities, —37.8 km s~^ for the 1989—1995 observations and —40.0 
km s~^ for the 2005—2010 observations, with no transition between them. 

This does not mean, however, that we have no idea of period of the suspected binary. 
Given the observing intervals, it would not be possible for a binary to go through a complete 
cycle of variability in a time less than the total observing interval of 22 years. Thus we have 
a minimum value for the period of the suggested binary orbit. Continued monitoring will 
help to further constrain this by revealing a period or increasing the minimum value. 
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Thus on the basis of the change in the radial velocities of IRAS 22272+5435, we make 
a tentative, but we think probable, identification of a binary companion. If we make the 
assumption that the long-term velocity difference is due to binary motion, we can then 
carry out a radial velocity solution and determine limiting binary parameters. For this 
we use the minimum values of 22 years and 1.3 km s~^ for the values of the period and 
observed velocity semi-amplitude, the latter being half of the difference in the average 
velocity of the two intervals when the 128.4 day period is removed from the data. Assuming 
a circular orbit, MppN=0.62 (which appears to be a typical post-AG B core raass 



based on models for a sta r with an initial main sequence mass of 2-3 Mq (IBlocker 



1995 



Vassiliadis fc Wood 



19941 )). and a range of inclinations, one finds from the mass function 



a range of possible masses for the secondary: M2 = 0.10 M0 (2=90°), 0.12 M0 (2=60°), 
0.22 MfT^ (z=30°), and 0.51 Mq (z=15°). From a mid-infrared imaging and modeling study. 



Ueta et al.l fl2001l ) determined an inclination of the torus of 2=25°±3°. Assuming that the 



binary orbit is coplanar with the equatorial density enhancement, this leads to a companion 
mass of 0.27±0.04 Mq. This is a very reasonable value for a secondary companion, and the 
assumption that it is a main-sequence star would place it at a spectral class of ~M4 V (Cox 
2000). For a circular orbit, this leads to a separation of 8.0 AU = 1700 R©, a value not only 
well outside the present radius of the PPN of R~100 R©, but one that would have been well 
outside it when the star was at the tip of the AGB, with R<500 R©. We emphasize that 
these are only preliminary, minimum values, based on the hypothesis that this long-term 
velocity difference is due to binary motion. For example, if this is a binary with a period 
of 34 years and the same semi-amplitude as above, the calculated value for the mass would 
increase to 0.32 M©, and if the velocity semi-amplitude were 1.8 km s~^, the mass would 
increase to 0.50 M©. Calculations ha ve shown that a binary w ith such parameters can 



form a very narrow-waist bipolar PN ( ISoker &: Rappaport 



2OOOD, although in this case the 



visible HST image shows the nebula to not be this extreme in shape. These nebulae are 
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all sur rounded by a larger halo, representing the earlier AGB mass loss. 



Huggins et al 



( 120101 ) calculate the effects of different parameters {M2, i, a) on the shaping of an initially 
spherical AGB wind to produce the observed halo in a PPN. They run a model for a binary 
with parameters similar to the preliminary values that we found for IRAS 22272+5435, for 
which they predict an approx imately round halo. This is consistent with that observed for 



this PPNe (iSahai et al. 



20071 . classify the halo as centrosymmetric with arc-like features). 
(Given an estimated distance of 1.9 kpc and a luminosity ratio of RiSxlO^ based on the 



M4 V spectral classification, there is no present hope of observing the companion directly.) 



3.3. Selection Effects and Limits on Binarity from the Null Results 

For the other six PPNs, no systematic change in velocity has occurred between the 
1991 — 1995 and the 2007—2010 observations, nor is any seen when we include additional 
published velocities. Period analyses show only the pulsational velocity variation. However, 
one needs to consider selection effects. Since we observed the brightest PPNs, they are not 
highly obscured. Thus, if they have a bipolar structure with an enhanced equatorial density 
region, they are probably biased towards having their equatorial planes close to the plane 
of the sky, thereby reducing their obscuration. This would reduce the observed orbital 
velocity, given the reasonable assumption that the binary orbit is in the equatorial plane. 

In spite of the uncertainty in the inclination of the equator (which mid-infrared images 
and spatial- kinematical observations can help solve), these null results still set useful limits 
on the binary parameter space. These are shown graphically in Figure |2l assuming that 
MppN=0.62 M0, circular orbits, and a conservative detection limit of K=2.0 km s^^. This 
shows that to remain undetected, a companion of 0.40 Mq must have an orbital period of 
P>3.5 yr if i>15° and P>24.5 yr if «>30°; an undetected companion of 0.25 Mq must have 
an orbital period of P>1.1 yr if 2>15°, P>8 yr if ^>30°, and P>23 yr if z>45°. Another 
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way to look at it is that M2>0.65 Mq is excluded for P<10 yr at z>15° and M2>0.30 M© is 
excluded for P<13 yr at i>30°. Lower-mass or longer-period companions than these would 
escape detection. 



EDITOR: PLACE FIGURE [2] HERE. 



We do know something about the orientation of the bipolar structure in these six PPNs 
based on HST images and on 2-D modeling in several cases with resolved mid-infrared 
images. HST images of IRA S 18095+2704 and 19475+3119 sug gest that they are viewed 



at some intermediate angles (lUeta et al. 



2000 



Sahai et al. 



20071 ). Models based in mid-IR 



images and, for IRAS 07134+1005, spatially-resolved molecular-line spectroscopy, have 
determined the approximate inclinations of the polar axis with r espect to the plane of 



the sky for two of the objects. For 



and 



br IRAS 17436+5003, i ~ 90° flMeixner et al. 



RAS 07134+1005, z ~ 80° flMeixner et al. 



20021) or z ~ 10° ( 



20041). 



GledhiU fc Yates 



20031 ). although for the latter of these objects one sees very disparate results. For IRAS 



22223+4327, a comparison of mid-infrared and HST visible image s implies an inclination 



that is closer to edge-on (90°) thai i pole-on ( 



the inclination is more uncertain (IGledhill et al. 



Clu 



be fc GledhiU 



2001 



20061 ). For IRAS 19500-1709, 



Volk, personal communication 2008). 



Therefore at least some of these appear to be inclined out of the plane of the sky with i > 
30°, and there does not appear to be an important bias of the sample to the plane of the 
sky. In this case the projection effects are less severe and the constraints on the mass or 
period of any undetected binary companion are more significant. Thus it appears that our 
non-detection of binarity in all six of these PPNs cannot be explained away as simply due 
to a low binary inclination to the plane of the sky. 
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DISCUSSION AND CONCLUSIONS 



This radial velocity study of seven PPNs shows only one has velocity variations, in 
addition to pulsation, that can tentatively be attributed to a binary companion. The 
presence of pulsation does complicate the search for a binary companion and makes it more 
difficult, but it does not invalidate the null result for these other six. As a counter example, 
we found for the related post-AGB object 89 Her that we could detect a binary companion 
(K=3.3 km s^^, P=290 d) even with a pulsating central star (K=1.6 km s^^, P=66 d; in 
preparation). Thus such binaries could be detected but were not. We used these null results 
to set limits on the properties of any undetected companions. The one tentative binary has 
a long period (>22 years) and probably a normal stellar mass companion (M2>~0.27 Mq). 



One might initially be surprised by this low binary fraction in light of t 



le discovery 



of a 



2003 



arge number of post— AGB bin aries by Van Winckel and collaborators ( IVan Winckel 



2007 



Van Winckel et al. 



20061 ). They find thus far that 27 out of a sample of 51 



post-AGB stars are spectroscopic binaries. However, these 51 post-AGB are not an unbiased 



sample but were chosen because they 



previously known post-AGB binaries ( IDe Ruyter et al 



possess several o 



the observed characteristics of 



20061 ). These post-AGB binaries 



are a distinctly different set of objects than the PPNs. They show a broad infrared excess 
(broad SED), indicating both hot and cool dust, and have abundance anomalies thought 



to be due to chemical fractionation of refractory e 



non-refractory elements by the star (jVan Winckel 



ements onto dust, with re-accretion of 



20031 ). These properties are attributed 



to the presence of a circumbinary disk. Most and perhaps all of these objects are binaries, 
with P fti 100—2600 days and e = 0.0—0.6. We would have detected such binaries but 
did not. In these post-AGB binaries, it is the binary that is thought to be responsible 



for forming and stabilizing the disk ( IVan Winckel et al. 



20061 ). The orbital periods of 



the shorter of these are of such that the systems would not accommodate within them a 



- 14 - 



large AGB star. Thus it appears tha t it is their binary nature that leads to their special 



properties (IWaelkens fc Waters 



|2004( ) and brings them to our attention due to their infrared 
excesses. PPNs do not share these properties, but rather display a clearly double-peaked 
SED with a much larger infrared excess, indicating a detached shell and much larger mass 
loss. The PPNs, at least the carbon-rich ones, have abundances in agreement with AGB 
nucleosynthesis, and none have the abundance anomalies seen in the post-AGB binaries. 



Also, the PPNs show a visible nebula which most of the post - AGB binaries do not (the Red 

De Marcol ( 2009 ). we conclude that the 



Rectangle is an exception ). Thus, in contrast to 
binary post-AGB objects in general represent a class of objects that are unlikely to evolve 
into PNs and therefore do not bear directly on the question of the shaping of the nebulae. 
PPNs, on the other hand, give every indication of being the immediate precursors of PNs. 

We can make the comparison instead to the binary central stars of PN s. Photometric 
searches indicate that 1 0—20 % of all PNs have a close (P<8 d) companion (IMiszalski et al. 



2009a : 



De Marco 



20091 ). With their short periods, it seems likely that the binary PNs 



formed through common envelope evolution in which the two stars did not merge. Might 
the PPNs be binaries, but presently in the common envelope stage? Since the common 



envelope stage is calculated to 



3e ver y short, on the order of the pre-common envelope 



orbital period ( iRicker fc Tamm 



20081 ). this cannot be the case, for it would be far too 
improbable to find six of our seven in this very short-lived phase. We cannot make a 
comparison with the fraction of PNs with a period in the range of 0.1 to 30 years, since this 
is observationally u nknown . Radial velo city studies of PNs with a resolution of 3 km s~^ 



have been initiated ( IDe Marco 



20061 ). However, these are complicated by the broad lines in 



the central stars and their variable winds, and no definitive results have been obtained. 

What do we know about the shapes of the binary PNs and how do they compare 
with the shapes of these PPNs? Based on a sample of 30 of these binary PNs with good 
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images, it has been determined that ~30% have nebulae with clear bipolar morphologies, 



and it is suggested that thi s number might 



other factors are included (IMiszalski et al. 



be as high as ~60% if inclination effects and 



2009bl ). This result is highly suggestive that 



a common envelope evolution without merger will commonly produce a bipolar nebula. 
However, this does not imply the inverse, that bipolar nebulae have a binary central star. 



All of our seven PPNs 
classification scheme of 
classification scheme of 



lave a bright cen tral s tar and would be class ified as SOLE in the 



Ueta et al 



Sahai et al 



(12000 



and 



Siodmiak et al 



(I2OO8I ). In the more detailed 



(120071 ) ■ six of the seven are classified: four as elongated, 
one as bipolar, and one as multipolar (see Tabled]). This suggests they each have an axis 
of symmetry that might arise from a equatorial density enhancement, and in several of the 
cases this enhancement is seen in the mid-infrared images. 

One is still left with the question of why only one of our seven PPNs shows evidence 
of being a binary, given that the binary fraction of stars is so muc h higher. The careful 



study 



of a sample of 164 solar-type (F7— G9 IV- V, V, VI) stars by 



Duquennoy fc Mayor 



(119911 ) finds ~50 % to be binaries. But this included visual binaries and common proper 
motion pairs, and resulted in a mean period of 180 years. If we restrict our comparison to 
the spectroscopic binaries, the fraction drops to ~25 % (with orbits) or ~33 % (including 
those detected to vary in velocity as binaries but without determined orbits). These results 
are based on high-precision velocities {a < 0.3 km s^^) over an observing range of up to 
13 years (average 8.6 years) and are a better comparison with our radial velocity sample. 
Given their higher precision and the absence of the pulsational variations which complicate 
the study of our stars, our tentative detection of one in seven (14 %) to be binaries does 
not appear to be anomalously low. 

Might these PPNs be binaries but with periods longer than 25 yr? These might not be 
detected in this radial velocity study, but they could still affect the shaping of the nebula. 
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although their effect would be reduced with increasing orbital period and separation. Might 
they be binaries but with low-mass (<0.25 Mq) companions? Our above limits on binarity 
do not exclude such companions. If the companion is a brown dwarf or a super- Jupiter 
planet, then it would escape detection in our program. While these can have a significant 
effect on the mass loss in certain cases and produce elliptical nebulae, it is estimated that 



planets will significantly affect mass loss in only 4-10% of AGB stars ( jLivio fc Soker 



20021) 



The results of this present radial velocity study provide the first direct test of the 
binary hypothesis in shaping PPNs, the direct precursors of PNs. While they do not rule 
it out, it seems fair to say that they do not support the binary hypothesis. Although this 
study has not answered the question of whether the shaping of PPNs and PNs is ultimately 
due to a binary companion, it has set significant constraints on the properties of a binary 
companion during the PPN phase. The lack of detection of a companion in six of the cases 
probably implies that any such companion either has a period that is very long (>25 yr) or 
a mass that is very low (<0.25 Mq). The effects of these on the mass loss and its shaping 
are obviously less than in the case of a higher mass, shorter period companion. These 
constraints can help guide future attempts to model the formation of the circumstellar 
density asymmetries with a binary companion. And of course they might not be binaries, 
and the asymmetric mass loss would then be due to something else. We know that they 
were pulsating during the previous AGB phase, and perhaps pulsation coupled with some 



other mechanism such as cool star spots (ISoker 



2000l ) is the mechanism responsible for the 



shaping. These results also do not appear to support the hypothesis that the inte nsive mass 



loss at the end of the AGB (the "superwind" ) is driven by a binary companion (iDe Marco 



20091 ). Since in these seven PPNs it is apparent that the envelopes are detached and the 
shaping of the nebula has started, these results might suggest two ways to form the shapes 
seen in PNs: (a) through common envelope evolution, as evidenced by the close binary 
nuclei of some PNs, or (b) through a non-common envelope process, which is occurring in 
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these PPNs. This latter process might involve a distant and/or low-mass companion or be 
due to a single, pulsating central star. 

This radial velocity study is continuing so that we can extend the temporal baseline in 
the search for evidence of even longer period binaries and seek to confirm the one tentative 
case. We have also begun a radial velocity study of several edge-on bipolar PPNs, in 
which we make the reasonable assumption that the binary orbit would be in the plane of 
the equatorial density enhancement. In such a case, we would see the full orbital velocity 
variations without suffering from an inclination effect. While in these cases the star is 
completely obscured from view in visible light, it is seen in the near-infrared and thus 
amenable to near-infrared spectroscopy. These PPNs have bipolar lobes and an obscuring 
dust lane, implying very strong shaping of the outflow. Rotation will also be investigated 
by comparing these edge-on cases with ones that are more nearly pole-on; since the pole-on 
ones are expected to appear as slow rotators, this comparison can give evidence in the 
edge-on cases of possible rotational spin-up or merger by a companion. 

We thank R. McClure and C. Waelkens for sharing their unpublished observations 
with us, D. Crabtree for making some observations for us, and D. Westpfahl for making 
some initial observations which showed the variability of one of these objects. We thank 
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22272-1-5435. BJH acknowledges the support of a University Research Professorship 
and a sabbatical leave from Valparaiso University and the hospitality of the Dominion 
Astrophysical Observatory during the initial stages of this research project. BJH also 
acknowledges support from the National Science Foundation (AST 9018032, 0407087, 
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Table 1. List of Program Objects 



IRAS ID 


HD 


Var. Star 


Other Names 


V(niag) 


Sp.T. 


Morph. Class.'' 


07134+1005 


56126 


CY CMi 


LS VI +10 15 


8.2 


F5 I 


Ec*(0.41),h(e) 


17436+5003 


161796 


V841 Her 




7.1 


G3 lb 


Ec*(0.41) 


18095+2704 




V887 Her 




10.4 


F3 lb 


Bc,h 


19475+3112 


331319 




LS H +31 9 


9.4 


F3 I 


Mc*(0.43),ps(in,s),h 


19500-1709 


187885 


V5112 Sgr 




8.7 


F3 I 




22223+4327 




V448 Lac 


DO 41288 


9.7 


GO la 


Ec*(0.43),h(a) 


22272+5435 


235858 


V354 Lac 




9.0 


G5 la 


Ec*(0.55),h(a) 



^According to the classification system of ISahai at ahl (|2007l ). The main classifications are E - 
elongated, B - bipolar, M - multipolar; c indicates closed lobes, * that the star is visible, ps the 
presence of point symmetry, and h the presence of a halo. For more details, see their paper. 



Table 2. Summary of Radial Velocity Observations 



IRAS ID 


Number of Observations 


<V_R> (km s~^ 






Vls_r(CO,OH)'= 


¥^(00,011)"^ Comments 




1991-95° 


2007-10'= 


Total^ 


1991-95"= 


2007-10"= 


Total^ 


(km s^"'^) 


(km s^"'^) 


(km s^"'^) 


07134+1005 


21 


18 


141 


88.0(0.8) 


86.3(0.5) 


87.5(0.3) 


14 


73 


88 ••• 


17436+5003 


59 


45 


178 


-53.1(0.2) 


-52.6(0.4) 


-52.6(0.2) 


11 


-35 


—54 longer P? 


18095+2704 


47 


29 


77 


-29.4(0.3) 


-30.5(0.3) 


-29.8(0.2) 


10 


-5 


-25 ■•• 


19475+3112 


38 


29 


71 


2.1(0.4) 


1.9(0.5) 


1.7(0.3) 


13 


18 


■•• 


19500-1709 


35 


13 


58 


14.5(0.5) 


13.8(0.7) 


13.9(0.3) 


12 


25 


13 ••• 


22223+4327 


34 


36 


81 


-40.5(0.4) 


-41.9(0.3) 


-41.3(0.2) 


11 


-30 


-42 • • ■ 


22272+5435 


34 


36 


155 


-37.6(0.3) 


-40.6(0.3) 


-39.4(0.2) 


12 


-28 


—40 longer P 



^The values in parentheses represent the uncertainties in the mean values. 
''Range in radial velocities, based on the total data set. 

"=00 or OH m o lecula r -line velocities in the loca l stand ard of rest (LSR) system. References: CO 



Likkel et al 



Bujarrabal et al 



(I1992II : iHrivnak fc Bieginel (|200,5h : OH 



Eder et al. 



(19M) 



|l99lh : 



Omont et al 



(|l993h : 



"^¥(00, OH) transformed from the LSR to the heliocentric system. 
''Based on our observations only. 



^Including additional data available from the literature or by personal communication as follows: IRAS 07134+1005 



Lebre et al 



( 1996 ) 



Barthes et al 



pOOOf ). Klochkova 



2003|); 17436+5003 



2009 p ersonal commu nication) 



Klochkova et al. 



(personal communication), 



R. McClure (personal communication), 




Van Winckel fc RevniersI (|2000l) . 



Klochkova et al 



Hrivnak fc Reddv 



( 2007h: 19500-1709 - C. Waelkens 



Van Winckel fc RevniersI pOOOh : 22272+5435 

iaood). 



to 
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Fig. la. — (Left) Radial velocity curves from our 1991 — 1995 observations, showing sine- 
curve fits to the periodic pulsational variations. (Right) Radial velocity curves showing the 
long-term velocity variations, including both our observations (filled circles) and those of 
others (open circles). Sample average error bars (±1 a) are shown for each PPN in the lower 
right-hand corner of the left panels. 
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Fig. lb. — Radial velocity curves, with symbols similar to Figure la. The dashed lines in the 
right-hand panel for IRAS 22272+545 represent the average velocities in the two observing 
intervals. 
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Fig. 2. — Plot showing limits on M2 and P for various inclinations (marked on the curves), 
assuming MppN=0.62 Mq and K<2.0 km s~^, based on the null results for binary detection 
for the other six PPNs. To remain undetected, any binaries present in the sample must have 
companion masses lower than or periods longer than the values of the curves at the various 
inclinations (i.e., they must lie below the curves). 



